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ABSTRACT 

We present results of a 1.1 mm deep survey of the AKARI Deep Field South (ADF-S) 
with AzTEC mounted on the Atacama Submillimetre Telescope Experiment (ASTE). 
We obtained a map of 0.25 deg 2 area with an rms noise level of 0.32-0.71 mJy. This 
is one of the deepest and widest maps thus far at millimetre and submillimetre wave- 
lengths. We uncovered 198 sources with a significance of 3. 5-15. 6a, providing the 
largest catalog of 1.1 mm sources in a contiguous region. Most of the sources are not 
detected in the far- infrared bands of the AKARI satellite, suggesting that they are 
mostly at z ^ 1.5 given the detection limits. We constructed differential and cumu- 
lative number counts in the ADF-S, the Subaru/XMM Newton Deep Field (SXDF), 
and the SSA 22 field surveyed by AzTEC/ASTE, which provide currently the tight- 
est constraints on the faint end. The integration of the best-fit number counts in the 
ADF-S find that the contribution of 1.1 mm sources with fluxes ^1 mJy to the cosmic 
infrared background (CIB) at 1.1 mm is 12-16%, suggesting that the large fraction 
of the CIB originates from faint sources of which the number counts are not yet con- 
strained. We estimate the cosmic star- formation rate density contributed by 1.1 mm 
sources with ^1 mJy using the best-fit number counts in the ADF-S and find that it is 
lower by about a factor of 5-10 compared to those derived from UV/optically-selected 
galaxies at z ~ 2-3. The fraction of stellar mass of the present-day universe produced 
by 1.1 mm sources with ^1 mJy at z ^ 1 is ~20%, calculated by the time integra- 
tion of the star-formation rate density. If we consider the recycled fraction of >0.4, 
which is the fraction of materials forming stars returned to the interstellar medium, 
the fraction of stellar mass produced by 1.1 mm sources decrease to <10%. 

Key words: submillimetre - galaxies: starburst - galaxies: evolution - galaxies: 
formation - galaxies: high-redshift 
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Over the past decade, millimetre and submillimetre ob- 
servations have shown that (sub) millimetre-bright galax- 
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ies (hereafter SMGs) hold important clues to galaxy evo- 
lution and the cosmic star formation history (|Blain et al.l 
2002, for a review). SMGs are highly obscured by dust, 
and the resulting thermal dust emission dominates the bolo- 
metric luminosity. The source of heating energy is domi- 
nated by vigorous star formation with star formation rates 
(SFRs) of several 100-1000 M© yr" 1 . Optical/near-infrared 
spectroscopy of a sample of SMGs with radio counter- 
parts revealed a median redshift of z ~ 2 for t he popula- 
tion dSwinbank et al.ll2004l : [ciiaDman et al.ll2005l). Recently 
SMGs at z > 4 have been confirmed (ICapak et al. [2008; 
ICoppin et ai1l2009l : iDaddi et all 120091 : iKnudsen et alJl2010h 
and ther e is now a spectros copically confirmed source at 
z = 5.3 (jRiechers et alJl2010h . 

Coupled w ith reports of high dynamical mass and 
gas mass (e.g., iGreve et a l. 2005; iTacconi et £11120061 ). it 
is suggested that SMGs are the progenitors of massive 
spher oidal galaxies observed during the ir formation phase 
(e.g., iLillv et al.l[l996l : ISmail et alJbooA . 

Mountin g evidence shows that the cosmic infrar ed back- 
ground (CIB IPuget et al.lll996l : iFixsen et al.lll998l ) at mil- 
limetre and submillimetre wavelengths is largely contributed 
by high-redshift galaxies (|Lagache et al.1 l2QQ5h . The CIB 
is the integral of unresolved emission from extragalactic 
sources and contains information on the evolutionary history 
of galaxies. While 850 /am surveys have resolve d 20-40% of 
the CIB into point sour c es in blank fields (e .g., lEales et"al1 
l200d : iBorvs et al.1 120031 : ICoppin et all 120061) and 50-100% 
in lensing cluster fields (e .g, iBlain et al. " 120021 : ICowie et al.l 



2002; Knudsen et al.l 120081), 1 mm blank field surveys have 
resolved only ~10% (e.g.. IGreve et al. _ 20041: [Laurent et aT 



2005; Malone v etaD 120051 : IScott et al.1 120081 : IPerera et al 
2008: IScott et alJl20ld ). A large portion of the CIB at mil- 



limetre and submillimetre wavelengths likely arises from 
galaxies with fainter flux densities. 

In conjunction with constraints from the CIB, the num- 
ber counts of SMGs are sensitive to the history of galaxy 
evolution at high redshift s. This requires constraining both 
the faint and bright end of the number counts, which in 
turn requires a suitable combination of small, deep sur- 
veys along with shallower wide-area surveys. Blank field 
surveys at millimetre and submillimetre wavelengths have 
been carried out with large bolometer arrays such as the 
Submillimetre Co mmon User Bolometer Array (SCUBA; 
iHolland et all 1 19991 ) on the James Clerk Maxwell Tele- 
scope (JCMT) (e.g ., ISmail et al.lll997l : iHughes et al.lfl998l : 
ICoppin et"aDl2006), the Max-Plank Millimetre Bolometer 
Array fMAMBO: iKrevsa et all Il998l) on the IRAM 30-m 
telescope and the Bolocam (|Glenn et al.lfl998l ) on the Cal- 
tech Subm il limetre Observatory (CSO) (e.g.. IGreve et al.l 
120041 . I2Q08I : iBertoldi et all 120071 : lLaurent et al.1 120051). the 
Larg e Apex BOlometer CAmera fLABOCA: Isiringo et aD 
2009) on the Atacama Pathfinder Experim ent (APEX) 
(e.g.. IWeiB et al.ll2009l : ISwinbank etaIll201Ch. and AzTEC 
(I Wilson et al.1 l2008al) on the JCMT (IScott et~ai1 120081 : 
IPerera et al.ll2008l : lAustermann et alJl2009l . l201Ch . However, 
the total area covered in existing surveys is still small 
(<1 deg 2 ) compared to the cosmic large scale structure, and 
substantial field-to-field variations can be seen in the pub- 
lished number counts. In addition, because of the limited 
depth of these surveys, the number counts of SMGs at faint 
flux densities (~1 mJy) are still not well constrained. 



We performed extensive surveys at 1.1 mm with AzTEC 
mounted on t he Atacama Submi ll imetr e Telescope Experi- 
ment (ASTE; lEzawa et al l 1 120041. 120081) in 2007 and 200 8 
(jWilson et al.ll2008bl : iTamura et al.ll2009l : IScott et alJboiOl ) . 
The ASTE is a 10-m submillimetre telescope located at 
Pampa la Bola in the Atacama desert in Chile. Some of 
the AzTEC/AS TE sources are followed up by submm/mm 
interferometers ([Hatsukade et al1l2010l : Ukarashi et al.l l2010. 
Tamura et al. in preparation). In this paper, we report 
on a deep blank field survey of the AKARI Deep Field- 
South (ADF-S). The ADF-S is a multi- wavelength deep 
survey field near the South Ecliptic Pole. It is known to 
be one of the lowest-cirrus region in the whole sky (100 
/im flux density of <0.5 MJy sr -1 and HI column den- 
sity of -5 x 10 19 cm" 2 ; ISchlegel et aD Il998l ). providing a 
window to the hi gh-redshift dusty unive rse. AKARI, an 
infrared satellite (Ma tsuhara et aO 120051 ). h as conducted 
deep surveys with the InfraRed Camera (IRC: lOnaka et all 
120071 ) at 2.4, 3.2, 4.1, 7 11, 15, 24 /im, an d with the Far- 
Infrared Surveyor fFIS: lKawada et al.ll2007l) at 65, 90, 140 
and 160 /im, down to the confusion limit (|Shirahata et al.l 
120091 ) . Multi- wavelength follow-up observations from the UV 
to the radio are under way. The full dat a sets of IR to 
subm illimetre bands with AKARI, Spitzer ([Clements et al.l 
2010), the Balloon-borne Large- Aperture Submillimetre 
Telescope (BLAST), Herschel Space Observatory, and the 
AzTEC/ASTE offers a unique opportunity to study the 
dusty galaxy population that contributes to the cosmic back- 
ground at IR-mm wavelengths. 



This paper presents the 1.1-mm map and source cata- 
log of the ADF-S. Together with the results from the Sub- 
aru/XMM Newton Deep Field ( SXDF; Ikarashi et a l. in 
preparation) and the SSA 22 field (jTamura et al.ll2009l ) sur- 
veyed by AzTEC/ASTE, we present statistical properties of 
the SMG population. Currently, these are the deepest wide- 
area survey ever made at millimetre wave lengths along wit h 
the AzTEC/ASTE GOODS-S survey (jScott et al.1 12OIOI ). 
providing the tightest constraints on number counts toward 
the faint flux density, albeit with lower resolution than has 
been typically employed to date. Comparisons with multi- 
wavelength data and statistical studies such as clustering 
analysis of this dataset will be presented in future papers. 



The arrangement of this paper is as follows. Sec- 
tion [2] summarises the observations of the ADF-S with 
AzTEC/ASTE. Section [3] outlines the data reduction and 
calibration details. In Section [4] we present the 1.1 mm map 
and the source catalogue. In Section [5] we derive number 
counts of the ADF-S, the SXDF, and the SSA 22 field, and 
compare them with other 1-mm wavelength surveys and lu- 
minosity evolution models. In Section[6]we estimate the con- 
tribution of 1.1 mm sources to the CIB. In Section [71 we con- 
strain the redshifts of the AzTEC sources using flux ratios 
between 1.1 mm and 90 /xm obtained with AKARI /FIS. In 
Section[8]we discuss the cosmic star formation history traced 
by 1.1 mm sources. A summary is presented in Section [9] 



Throughout the paper, we adopt a cosmology with 
Ho 70 km s" 1 Mpc" 1 , h 0.7, Q M 0.3, and Q A 0.7. 
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Figure 1. Signal map of the ADF-S. White contours represent the 30% (outer contour) and 50% (inner contour) coverage regions 
(corresponding to rms noise levels of 0.55 and 0.71 mJy beam -1 ), respectively. The 198 sources detected with ^3.5cr in the 30% coverage 
region are indicated by green circles with a 30" diameter, and numbered in order of significance (Table [2]). The pixel size is 3" x 3". 



2 OBSERVATIONS 

The center region of the ADF-S was observed with AzTEC 
on the ASTE. The observations were made from September 
16 to October 14 in 2007, and from August 4 to December 
21 in 2008. The operations were carried out remotely from 
the ASTE operation rooms in San Pedro de Atacama, Chile, 
and in Mitaka, Japan through the network observation sys- 
tem N-COSMOS3 developed b y the National Astrono mical 
Observatory of Japan (NAOJ) (jKamazaki et al.ll2005h . The 
full width at half maximum (FWHM) of the AzTEC detec- 
tors on the ASTE is ~30" at 1.1 mm (270 GHz), and the 
field of view of the array is roughly circular with a diame- 
ter of 8'. During the observing run, 107 and 117 out of the 
144 AzTEC detectors were operational in 2007 and 2008, 
respectively. 

In order to maximize the observing efficiency, we used 



a lissajous scan pattern to map the field. We chose a high 
maximum velocity of 300" s _1 to mitigate low- frequency at- 
mospheric fluctuations. The lissajous scan pattern provides 
~20 / x 20' coverage on the sky, in which the central ^12' x 12' 
area is nearly uniform, with decreasing integration time from 
the inside to the outside. We covered the ADF-S with seven 
different field centers to make the noise level of the entire 
map uniform. We obtained a total of 319 individual observa- 
tions under conditions where the atmospheric zenith opacity 
at 220 GHz as monitored with a radiometer at the ASTE 
telescope site was T220 GHz = 0.02-0.1. The total time spent 
on- field was ^216 hours. Uranus or Neptune was observed 
at least once a night in raster-scan mode to measure each 
detector's point spread function (PSF) and relative position, 
and to determin e the flux conversion factor (FCF) for abso- 
lute calibration ([Wilson et al. 2008a). Pointing observations 
with the quasar J0455-462 that lies ^7 degree from the field 
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center were performed every two hours, bracketing science 
observations. A pointing model for each observing run in 
2007 and in 2008 is constructed from these data, and we 
make corrections to the telescope as trometry. The resu lting 
pointing accuracy is better than 3" (Wil son et al.ll2008bl ). 



3 DATA REDUCTION 

The data were reduced in a manner identical to that de- 
scribed in IScott et"al ] <|2Q10h . We used a principle compo- 
nent analysis (PCA) to remove the low-fre quency atmo- 
spheric signal from th e time-stream data ( Laurent et al.l 
120051 : IScott et alJl20Q8h . The PCA method AC couples the 
bolometer time-stream data, making the entire map and 
point-source kernel have mean of zero. The cleaned time- 
series data is projected into map space using 3" x 3" pixels 
in R.A.-Dec, and the individual observations are co-added 
into a single map by weighted averaging. We also create 100 
noise reali zations by 'ja ckknifing' the time-series data as de- 
scribed in IScott et al.1 (|2008h . These noise maps represent 
realizations of the underlying random noise in the map in 
the absence of sources (both bright and confused) and are 
used throughout this paper to characterize the properties of 
the map and source catalogue. The point source profile is 
affected by the PCA method since the faint point sources 
with low spatial frequencies are also attenuated. The PCA 
method makes the mean of the map zero, causing negative 
side lobes around the peak of the point source profile. We 
trace the effects of PCA and other process in the analysis 
on the point source profile; this 'point source kernel' is used 
to optimally filter the coadded map and the 100 noise real- 
izations for the detection of point sources. A 2-dimensional 
Gaussian fitting to the point source kernel gives a FWHM 
of35'/9. 



Table 1. Map properties in the 50% and the 30-50% coverage 

rorn r\n c 



regions 



Coverage 


Area 


Noise level 


Source^ 


False Detections^ 


(%) 


(arcmin 2 ) 


(mJy beam -1 ) 




50% 


709 


0.32-0.55 


169 


4.9 ±0.22 


30-50% 


200 


0.55-0.71 


29 


1.4 ±0.12 



a Number of sources (^3.5cr). 

^ Number of false detection (^3.5<j). 



RMS Noise Level (mJy beam'1) 

0.6 0.8 1.0 
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Figure 2. Noise map of the AzTEC/ASTE ADF-S. The outer 
and inner contours indicate the 30% and 50% coverage region, 
respectively. 



4 MAP AND SOURCE CATALOGUE 
4.1 Map 

The signal map and the corresponding noise map are shown 
in Figures [1] and [21 respectively. The dashed curves represent 
regions with ^30% (outer contour) and ^50% (inner con- 
tour) of the maximum weights (hereafter called as 30% and 
50% coverage region). The area and noise level are 709 and 
200 arcmin 2 , and 0.32-0.55, 0.55-0.71 mJy in the 50% and 
30-50% coverage region, respectively (Table [T]). This survey 
is confusion- limited, where the 5cr-confusion limit estimated 
bv lTakeuchi &; Ishi i (2004) is 4.4 mJy using the point source 
kernel and the differential number counts in the ADF-S de- 
rived in § 15.11 

Fig. [3] shows the distribution of flux values in the map, 
compared to that averaged over the 100 noise realizations 
within the 50% coverage region. The result of a Gaussian fit 
to the averaged noise map is superimposed in Fig. [3j The 
presence of real sources in the map makes excess of both 
positive and negative valued pixels over the histogram of 
the noise map, since the signal map is created to have a 
mean of zero. This fit deviates from the distribution of pixel 
values at high positive and negative fluxes because the map 
is not uniform over the entire region, with the outer region 
being slightly noisier. 





Signal Map 




\ Noise Map 


: pT / 


Gaussian Fit 




, , , u : Ah 



Flux Density [mJy beam" 1 ] 



Figure 3. Distribution of flux density in the real map (solid 
histogram) and averaged over the 100 noise realizations (dashed 
histogram) in the 50% coverage region. The result of a Gaussian 
fit to the flux distribution of the noise maps is plotted as a dashed 
curve. 
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Figure 4. Number of false detections expected above given 
signal-to-noise ratio thresholds calculated in the 50% and 30-50% 
coverage regions. The error bars are lcr Poisson uncertainties. 

4.2 Source Catalogue 

Source extraction is performed on the SN map using a crite- 
rion of ^3.5cr. The source positions are determined by flux- 
squared weighting on pixels within a 15" radius of the nom- 
inal peak. We detect 198 and 169 sources with a significance 
of 3.5-15.6cr in the 30% and 50% coverage regions, respec- 
tively. The source catalog is given in Table [2] where both the 
observed and the deboosted flux densities (§ I4.6J) are listed. 
169 sources (ADFS-AzTECl-169) detected within the 50% 
coverage region, the deeper and more uniform coverage re- 
gion, are listed first, followed by the remaining 29 sources 
(ADFS-AzTEC170-198) detected outside the 50% coverage 
region. 



Figure 5. Survey completeness for the 50% coverage region. The 
error bars are lcr from the binomial distribution. 
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4.3 False Detections 

Monte Carlo simulations are carried out to estimate the 
number of spurious sources due to positive noise fluctua- 
tions. We conduct the standard source extraction on the 
100 synthesised noise realizations, and count the number of 
'sources' above given S/N thresholds in steps of 0.5cr. Fig. [4] 
shows the average number of false detections as a function of 
S/N. The expected number of false detections in our ^3.5cr 
source catalog is ^4-5 and ~l-2 in the 50% and 30%-50% 
coverage region, respectively. 



4.4 Completeness 

The survey completeness is computed by injecting simulated 
point sources of known flux densities into the real signal map 
one at a time. The input positions are randomly selected 
within the 50% coverage region, but are required to be out- 
side a 20" radius from a real source in the map to avoid 
blending. When a simulated source is extracted within 20" 
of its input position with S/N ^ 3.5, the source is considered 
to be recovered. We repeat this 1000 times for each flux bin, 
and compute the fraction of output sources to input sources. 
The completeness as a function of intrinsic flux density is 
shown in Fig. [5] The error bars are the 68% confidence in- 



Figure 6. Probability that a source is detected outside an angu- 
lar distance from its true position. The probability is calculated 
for sources with 3.5 < SN < 4.0, 4.5 < SN < 5.0, and 5.5 < 
SN < 6.0. The horizontal dashed lines indicate 68.3% and 99.5% 
confidence interval s. The shaded regions represent theoretical pre- 
dictions derived in llvison et al.l (|2007l ) for the SN ranges. 

tervals from the binomial distribution. The completeness is 
about 50% at a flux density of 2.0 mJy. 

4.5 Positional Uncertainty 

The positional error is an important indicator when iden- 
tifying other- wavelength counterparts. The positional un- 
certainties for detected sources are calculated in a similar 
manner for the completeness calculation. Fake sources of 
known flux densities (ranging from 1 to 10 mJy in steps of 
0.25 mJy) are injected into the real signal map one at a time. 
The input positions are selected randomly outside a 20" ra- 
dius from the real sources. Source extractions with S/N ^ 
3.5 are performed within 20" search radius from the input 
positions. This procedures are repeated 1000 times for each 
flux bin, and measure the angular distances between the in- 
put and output positions. The probability that a source is 
extracted outside an angular distance from its input posi- 
tion is shown in Fig. [6] for different S/N ranges. The 68.3% 
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Table 2. AzTEC/ASTE ADF-S 3.5cr Source Catalog. AzTECl-169 are detected in the 50% coverage region and AzTEC170-198 are 
detected in the region with coverage 30-50%. The columns give (1) Source name; (2) ID; (3) Right ascension; (4) Declination; (5) 
Observed flux density and la error; (6) Deboosted flux density and 68% confidence level; (7) Signal to noise ratio. 
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J044452. 


,99- 


-533810. 


.5 


ADFS- 


-AzTECll 


4 


44 


52, 


,99 


-53 


38 


10, 


.5 


4, 


,8 ± 0. 


,4 


4, 


. + 0.5 
' -0.5 


11.4 


AzTEC 


J044733. 


,99- 


-531837. 


.5 


ADFS- 


-AzTEC12 


4 


47 


33 


,99 


-53 


18 


37. 


.5 


4. 


,0 ± 0. 


.4 


3. 


O + 0.4 
— 0.4 


11.0 


AzTEC 


J044736. 


.90- 


-532527. 


.5 


ADFS- 


-AzTEC13 


4 


47 


36 


,90 


-53 


25 


27. 


.5 


4. 


,5 ± 0. 


,4 


4. 


O + 0.5 
' — 0.5 


10.7 


AzTEC 


J044442. 


,42- 


-534122, 


,9 


ADFS- 


-AzTEC14 


4 


44 


42, 


,42 


-53 


41 


22. 


,9 


4, 


,8 ± 0, 


,4 


4, 


k+0.5 
' — 0.5 


10.6 


AzTEC 


J044844. 


.20- 


-531443. 


,4 


ADFS- 


-AzTEC15 


4 


48 


44 


,20 


-53 


14 


43. 


,4 


5. 


,0 ± 0. 


,5 


4. 


7+0.5 
' -0.5 


10.5 


AzTEC 


J044802. 


.69- 


-531712. 


,7 


ADFS- 


-AzTEC16 


4 


48 


02 


,69 


-53 


17 


12. 


,7 


4. 


,2 ± 0. 


,4 


3. 


q + 0.4 
' — 0.5 


10.5 


AzTEC 


J044743. 


,53- 


-531541, 


,8 


ADFS- 


-AzTEC17 


4 


47 


43. 


,53 


-53 


15 


41. 


,8 


4, 


,2 ± 0. 


,4 


4, 


o + °-? 

-0.5 


10.3 


AzTEC 


J044753. 


,52- 


-532331. 


,4 


ADFS- 


-AzTEC18 


4 


47 


53 


,52 


-53 


23 


31. 


,4 


4. 


,1 ± 0. 


.4 


3. 


Q + 0.5 
— 0.4 


10.3 


AzTEC 


J044543. 


,44- 


-532524. 


.3 


ADFS- 


-AzTEC19 


4 


45 


43 


,44 


-53 


25 


24. 


.3 


3. 


,8 ± 0. 


,4 


3. 


u -0.4 


10.1 


AzTEC 


J044813. 


,58- 


-531859. 


.6 


ADFS- 


-AzTEC20 


4 


48 


13, 


,58 


-53 


18 


59, 


,6 


4, 


,1 ± 0. 


,4 


3, 


9 + 0-5 
' — 0.5 


10.0 


AzTEC 


J044658. 


.41- 


-531534. 


,1 


ADFS- 


-AzTEC21 


4 


46 


58. 


,41 


-53 


15 


34. 


,1 


5. 


,0 ± 0. 


,5 


4. 


7+0.6 
' ' —0.5 


9.9 


AzTEC 


J044455. 


,57- 


-533423. 


.0 


ADFS- 


-AzTEC22 


4 


44 


55 


,57 


-53 


34 


23. 


.0 


3. 


.9 ± 0. 


,4 


3. 


7+0.5 
' ' -0.4 


9.8 


AzTEC 


J044540. 


,95- 


-533318. 


,2 


ADFS- 


-AzTEC23 


4 


45 


40. 


,95 


-53 


33 


18, 


,2 


3, 


,2 ± 0, 


,3 


3, 


w — 0.4 


9.7 


AzTEC 


J044522. 


.00- 


-532700. 


.9 


ADFS- 


-AzTEC24 


4 


45 


22 


,00 


-53 


27 


00. 


.9 


3. 


.6 ± 0, 


.4 


3. 


— 0.4 


9.7 


AzTEC 


J044659. 


.94- 


-531907. 


,3 


ADFS- 


-AzTEC25 


4 


46 


59. 


,94 


-53 


19 


07, 


,3 


4, 


,0 ± 0. 


,4 


3, 


7+0.5 
' ' -0.5 


9.7 


AzTEC 


J044629. 


.83- 


-533345. 


.0 


ADFS- 


-AzTEC26 


4 


46 


29, 


,83 


-53 


33 


45, 


.0 


3, 


,3 ± 0. 


,4 


3, 


1+0.4 
— 0.4 


9.3 


AzTEC 


J044658. 


.48- 


-532327. 


.8 


ADFS- 


-AzTEC27 


4 


46 


58. 


,48 


-53 


23 


27. 


.8 


3. 


.6 ± 0. 


.4 


3. 


— 0.5 


9.1 


AzTEC 


J044612. 


,33- 


-532743. 


,6 


ADFS- 


-AzTEC28 


4 


46 


12. 


,33 


-53 


27 


43, 


,6 


2, 


,9 ± 0, 


,3 


2, 


7+0.4 
' ' -0.4 


9.1 


AzTEC 


J044714. 


,90- 


-531738. 


,9 


ADFS- 


-AzTEC29 


4 


47 


14. 


,90 


-53 


17 


38, 


,9 


3, 


,7 ± 0. 


,4 


3, 


•4 + °1 

— 0.5 


9.1 


AzTEC 


J044805. 


.34- 


-532433. 


.2 


ADFS- 


-AzTEC30 


4 


48 


05. 


,34 


-53 


24 


33. 


.2 


4. 


,4 ± 0. 


,5 


4. 


1+0.6 
— 0.5 


9.0 


AzTEC 


J044351. 


.90- 


-534443. 


,3 


ADFS- 


-AzTEC31 


4 


43 


51. 


,90 


-53 


44 


43, 


,3 


4, 


,8 ± 0. 


,5 


4, 


c + 0.6 
' -0.6 


8.8 


AzTEC 


J044553. 


,57- 


-533520. 


.6 


ADFS- 


-AzTEC32 


4 


45 


53, 


,57 


-53 


35 


20, 


.6 


2. 


,8 ± 0, 


,3 


2, 


^ — 0.4 


8.7 


AzTEC 


J044706. 


.39- 


-531612. 


.5 


ADFS- 


-AzTEC33 


4 


47 


06. 


,39 


-53 


16 


12. 


.5 


4. 


,0 ± 0. 


,5 


3. 


7+0.5 
' -0.5 


8.7 


AzTEC 


J044538. 


,56- 


-532827. 


,6 


ADFS- 


-AzTEC34 


4 


45 


38. 


,56 


-53 


28 


27, 


,6 


2, 


,8 ± 0. 


,3 


2, 


.6+°1 

-0.4 


8.4 


AzTEC 


J044853. 


,11- 


-531557. 


,9 


ADFS- 


-AzTEC35 


4 


48 


53. 


,11 


-53 


15 


57, 


,9 


4. 


,1 ± 0. 


,5 


3, 


7+0.5 
' ' -0.6 


8.2 


AzTEC 


J044556. 


.43- 


-533929. 


,7 


ADFS- 


-AzTEC36 


4 


45 


56. 


,43 


-53 


39 


29. 


,7 


3. 


,5 ± 0. 


,4 


3. 




8.1 


AzTEC 


J044348. 


,44- 


-534314, 


,2 


ADFS- 


-AzTEC37 


4 


43 


48. 


,44 


-53 


43 


14. 


,2 


4. 


,3 ± 0. 


,5 


3, 


9+0.6 

'^ — 0.6 


8.0 


AzTEC 


J044824. 


,87- 


-531501. 


.3 


ADFS- 


-AzTEC38 


4 


48 


24, 


,87 


-53 


15 


01, 


.3 


3, 


,5 ± 0. 


,4 


3, 


— 0.5 


8.0 


AzTEC 


J044609. 


.02- 


-532710. 


.8 


ADFS- 


-AzTEC39 


4 


46 


09 


,02 


-53 


27 


10. 


.8 


2. 


.6 ± 0. 


,3 


2. 


— 0.4 


8.0 


AzTEC 


J044534. 


,14- 


-533939. 


,6 


ADFS- 


-AzTEC40 


4 


45 


34. 


,14 


-53 


39 


39, 


,6 


3, 


,4 ± 0. 


,4 


3, 


1+0.5 
— 0.5 


7.9 


AzTEC 


J044354. 


,85- 


-533928. 


,8 


ADFS- 


-AzTEC41 


4 


43 


54. 


,85 


-53 


39 


28, 


,8 


3, 


,9 ± 0, 


,5 


3, 


-0.5 


7.9 


AzTEC 


J044630. 


.39- 


-533154. 


,7 


ADFS- 


-AzTEC42 


4 


46 


30 


,39 


-53 


31 


54. 


,7 


2. 


,7 ± 0. 


,3 


2. 


. + 0.4 
' — 0.4 


7.8 


AzTEC 


J044607. 


,77- 


-532813. 


.4 


ADFS- 


-AzTEC43 


4 


46 


07, 


,77 


-53 


28 


13, 


.4 


2, 


,5 ± 0, 


,3 


2, 


q + 0.4 
,J -0.3 


7.8 


AzTEC 


J044643. 


,79- 


-533002. 


.4 


ADFS- 


-AzTEC44 


4 


46 




, ( y 


-53 


on 
oU 


no 
uz. 


.4 


3, 


,0 ± 0. 


,4 


2. 


7+0.4 
' ' —0.5 


7.6 


AzTEC 


J044451. 


,57- 


-533544. 


.1 


ADFS- 


-AzTEC45 


4 


44 


51 


,57 


-53 


35 


44. 


.1 


3. 


,1±0, 


,4 


2. 


9 + 0.5 


7.5 


AzTEC 


J044420. 


,88- 


-534011. 


,4 


ADFS- 


-AzTEC46 


4 


44 


20, 


,88 


-53 


40 


11, 


,4 


3, 


,4±0, 


,5 


3, 


1+0.5 

-"--o^ 


7.3 


AzTEC 


J044541. 


,74- 


-533445. 


,2 


ADFS- 


-AzTEC47 


4 


45 


41. 


,74 


-53 


34 


45, 


,2 


2, 


,3±0, 


,3 


2, 


1+0.4 
• x -0.4 


7.3 


AzTEC 


J044437. 


,74- 


-534637. 


.8 


ADFS- 


-AzTEC48 


4 


44 


37. 


,74 


-53 


46 


37. 


.8 


3. 


.8±0. 


,5 


3. 


4 +o.6 

'^-0.6 


7.2 


AzTEC 


J044528. 


,39- 


-533715. 


.8 


ADFS- 


-AzTEC49 


4 


45 


28. 


,39 


-53 


37 


15, 


.8 


2, 


,4±0. 


,3 


2, 


9 + 0.4 


7.1 


AzTEC 


J044456. 


,63- 


-533616. 


,2 


ADFS- 


-AzTEC50 


4 


44 


56. 


,63 


-53 


36 


16, 


,2 


2, 


,8±0, 


,4 


2, 


. + 0.4 
,0 -0.5 


7.1 


AzTEC 


J044623. 


,49- 


-533552. 


.4 


ADFS- 


-AzTEC51 


4 


46 


23 


,49 


-53 


35 


52. 


.4 


2. 


,6±0, 


,4 


2. 


4+0.4 


6.9 


AzTEC 


J044344. 


,41- 


-534325. 


,9 


ADFS- 


-AzTEC52 


4 


43 


44. 


,41 


-53 


43 


25, 


,9 


3, 


,8±0, 


,6 


3, 


4+0.6 

-0.6 


6.9 


AzTEC 


J044742. 


,49- 


-531951. 


,2 


ADFS- 


-AzTEC53 


4 


47 


42. 


,49 


-53 


19 


51, 


,2 


2, 


,4±0, 


,4 


2, 


9 + 0.4 
z -0.4 


6.9 


AzTEC 


J044721. 


,15- 


-532659. 


.2 


ADFS- 


-AzTEC54 


4 


47 


21 


,15 


-53 


26 


59. 


.2 


2. 


.9±0. 


,4 


2. 


g + 0.5 
°-0.5 


6.8 


AzTEC 


J044423. 


,80- 


-533413. 


,5 


ADFS- 


-AzTEC55 


4 


44 


23. 


,80 


-53 


34 


13, 


,5 


3, 


,7±0, 


,6 


3, 


o + 0.6 
'^-0.6 


6.8 
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Table 2 - continued 





Name 








ID 




R.A. 




Dec. 






□bserved 


^deboosted 


S/N 


















■ 


m 


( c 


• / //j 






(mJy) 




( 


mJy) 




AzTEC 


J044435. 


35- 


OO^O^U. 


u 


Jr\ur 0- 


A7TFP^fi 


4 


44 


35. 


35 


—53 


43 


46. 


u 


Q 
O. 


,1 ± 0. 


.5 


2. 


g + 0.5 
°-0.5 


7 

u. / 


AzTEC 


J044611. 


,47- 




c 

.0 


Jr\ur 0- 


•l\Zi 1 LiVjO I 


4 


46 


11. 


,47 


—53 


21 


52, 


Q 
,0 


Q 
O. 


,5 ± 0. 


.5 


3. 


-1+0.5 
-0.6 


7 
u. / 


AzTEC 


T04451 6 


,36- 




Q 
, y 


ADFS- 


A^TFP^iR 

-Jr\Zi X H/^dO 


4 


45 


16. 


,36 


—53 


33 


09, 


Q 


2 


.2 ± 0. 


,3 


2. 


+0.4 
u -0.4 


7 
u. / 


AzTEC 


J044639. 


,12- 


uozuio. 


,u 


ADFS- 


A^TFP^iQ 

-Jr\Zi J_ HjkjOcj 


4 


46 


39. 


.12 


—53 


25 


18, 


,u 


2 


,5 ± 0. 


,4 


2, 


2 +o.4 

-0.4 


7 
u. / 


AzTEC 


J044527. 


,26- 




,u 


ADFS- 


AyTFPfiO 


4 


45 


27. 


,26 


—53 


35 


18, 


,u 


2 


,1 ± 0. 


,3 


1, 


g+0.4 
,y -0.4 


7 
u. / 


AzTEC 


J044858. 


,01- 


-^1 ^94 


3 


ADFS- 


A?TFPfi1 

■nZ X H/ v^U L 


4 


48 


58. 


.01 


—53 


15 


24 


3 


3 


,5 ± 0. 


,5 


3 


9 +0.6 
'^-0.6 


6.6 


AzTEC 


J044722. 


,55- 




Q 
,0 


ADFS- 




4 


47 


22. 


,55 


—53 


18 


25, 


c 

,0 


2 


,5 ± 0, 


,4 


2, 


3+0.5 

,o -0.4 


u.u 


AzTEC 


T0444nn 

o \j j. j. ±yj yj , 


,63- 


e;Q499» 


2 


ADFS- 




4 


44 


00. 


,63 


—53 


42 


28, 


2 


Q 




.3 ± 0, 


,5 


2, 


9+0-5 
,y -0.6 


u.u 


AZ 1 JZjv^ 


J044813. 


,94- 




7 


A F>T7Q 
±\LJr 0- 




A 
^± 


48 


13. 


.94 


—53 


22 


44 


7 


Q 
O 


.0 ± 0, 


,5 


2, 


7 +0.5 
• ' -0.5 


u.u 


AzTEC 


J044822. 


,25- 


c;qi 490 


n 

,u 


ADFS- 




4 


48 


22. 


.25 


—53 


14 


20, 


n 

,u 


2 


,9 ± 0, 


,4 


2, 


6 +0.5 
u -0.5 


U.u 


AzTEC 


J044503. 


,53- 




7 


Jr\ur 0- 




4 


45 


03. 


,53 


—53 


28 


16, 


7 


2 


,8 ± 0. 


,4 


2, 


,u -0.5 


U.u 


AzTEC 


J044826. 


.06- 


-531413. 


7 


ADFS- 




4 


48 


26. 


.06 


—53 


14 


13 


7 


2 


,9 ± 0. 


,4 


2, 


6 +0.5 
u -0.5 


6.5 


AzTEC 


J044857. 


,68- 


c;qi Q97 


n 

,u 


ADFS- 




4 


48 


57. 


,68 


—53 


18 


27, 


n 

,u 


Q 
O 


,5 ± 0. 


,5 


3 


1+O.6 
-0.6 


u.u 


AzTEC 




.20- 




2 


ADFS- 


A 7 TFPfiQ 


4 


45 


56. 


.20 


—53 


25 


41, 


2 


2 


.2 ± 0. 


4 


2, 


+0.4 


9 
u. z 


AZ 1 H/V^y 


J044603. 


,17- 


kqoo/17 
"000(54 I . 


, 




A 7 TFP7fl 

■AZ ± H/V^y < U 


/i 


46 


03. 


.17 


—53 


38 


47 


, 


9 
z. 


.6 ± 0, 


,4 


2, 


•°-0.4 


9 
u. z 


AzTEC 


J044358. 


,65- 


000 101, 


2 


ADFS- 


A7TFP71 
-l\Zi L Hjkj 1 L 


4 


43 


58. 


.65 


—53 


37 


31, 


2 


Q 
O. 


,1 ± 0, 


,5 


2, 


7 +0.5 
• ' -0.6 


9 
u. z 


AzTEC 


J044717. 


,48- 




c 
,0 


ADFS- 


A7TFP79 


4 


47 


17. 


.48 


—53 


26 


14, 


c 

,0 


2 


.5 ± 0, 


4 


2, 


9 +0.5 
,z -0.4 


u. 1 


az _L Hjkj 


J044428. 


,74- 




/I 


A F>T7Q 
±\LJr 0- 




/i 


44 


28. 


.74 


—53 


41 


29 


A 


9 
z. 


,8 ± 0. 


,5 


2, 


4+0.5 

'^-0.5 


u.u 


AzTEC 


J044747, 


.96- 


OOlOUu. 


, 


ADFS- 


A7TFP74 

■AZ 1 H/ 1 


4 


47 


47. 


,96 


—53 


13 


05, 


, O 


Q 
O 


,1 ± 0, 


,5 


2, 


7 +0.6 
' ' -0.6 


u.u 


AzTEC 


T044545 


,17- 


uoooiu. 


Q 
,0 


ADFS- 


A 7 TFP7 t ; 


4 


45 


45. 


,17 


—53 


38 


15, 


Q 
,0 


2 


.2 ± 0. 


4 


1 


g + 0.4 

,y -0.4 


u.u 


AZ 1 H/V^y 


J044502. 


,41- 


KQOQO/1 






A 7 TFP7fi 


A 

^± 


45 


02. 


.41 


—53 


39 


34 


A 

.^i 


9 
z. 


,5 ± 0, 


,4 


2, 


9 +0.5 
'^-0.5 


u.u 


AzTEC 


J044554. 


,73- 


c;Qqo9Q 
oooozy, 


,U 


ADFS- 


A7TFP77 


4 


45 


54. 


,73 


—53 


38 


29, 


,u 


2 


,3 ± 0. 


,4 


2, 


+0.5 


u.y 


AzTEC 


J044744. 


,88- 


-^1 R99 


,8 


ADFS- 


- AZ X H/ 1 O 


4 


47 


44. 


.88 


—53 


16 


22 


,8 


2 


,3 ± 0. 


,4 


2, 


+0.4 
U -0.5 


5.9 


az ± Hjkj 


J044614. 


,68- 


"OOZDOl, 


,0 




A 7 TFP7Q 


/i 


46 


14. 


.68 


—53 


26 


31 


,0 


9 
z. 


.0 ± 0, 


,3 


1 


7+0.4 
' ' -0.4 


u.y 


AzTEC 


J044635. 


,07- 


oozyou. 


7 


ADFS- 


A 7 TFPSn 

-AZ X H/V^OU 


4 


46 


35. 


,07 


—53 


29 


50, 


7 


2 


,1 ± 0, 


,4 


1 


g+0.4 


u.o 


AzTEC 


J044409. 


,74- 




4 


ADFS- 


-AZ X H/ V^O L 


4 


44 


09. 


.74 


—53 


46 


07, 


4 


3 


,1 ± 0. 


,5 


2, 


7 +0.6 
' ' -0.6 


5.8 


A 7 TPP 


J044432. 


,70- 




Q 


A F>T7Q 
±\LJr 0- 


A 7 TFPS9 
-AZ ± EjVjOZ 


/i 


44 


32. 


.70 


—53 


44 


19 


Q 


9 
z. 


,7 ± 0, 


,5 


2, 


4+0.6 

-0.5 


u.o 


AzTEC 


J044429. 


,65- 


uo^uuy . 


4 


ADFS- 


A 7 TFP8^ 

-AZ ± H/vyOO 


4 


44 


29. 


,65 


—53 


46 


59, 


4 


Q 
O 


,1 ± 0, 


,5 


2, 


7 +0.6 
' ' -0.6 


u.o 


A 7 TPP 


J044531. 


.48- 


KQoni q 
"Ooouiy, 


1 

, 1 


A F>T7Q 
±\LJr 0- 


-AZ ± H/V^o^ 


/i 


45 
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Table 2 - continued 
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Table 2 - continued 





Name 








ID 




R.A. 






Dec. 




s 
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^deboosted 


S/N 
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(mJy) 
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confidence interval for a 3.5cr source is ^9". We compare 
the positional uncertainty with theoreti cal predictions fo r 
uncorrelated Gaussian noise derived bv llvison et al.l (|2007l ). 
Shaded regions in Fig. [6] represent the theoretical probability 
distributions for sources with 3.5 < S/N < 4.0, 4.5 < S/N 
< 5.0, and 5.5 ^ S/N < 6.0 from top to bottom. The posi- 
tional uncertainties of AzTEC sources are broad compared 
to the theoretical predictions. It is possible that the ADF-S 
map is confused and the confusion noise affects source posi- 
tions. 



4.6 Flux Deboosting 

When dealing with a low S/N map, we need to consider the 
effect that flux densities of low S/N sources are boosted to 
above detection thresholds, due to the stee p slope of num- 
ber counts in the flux ran ge of our map ([Murdoch et al.l 
Il973l : iHogg fc Turner! [l998h . We correct for the flux boost- 
ing effect to estimate intrinsic flux densities of sources 
based on the Bayesian estimation (|Coppin et al. 2005. 120061 : 
lAustermann et al.|[2QQ9l . l2010h . 

To calculate the probability distribution of its intrinsic 
flux density (posterior flux distribution; PFD), we use the 
best-fit differential number counts in the ADF-S derived in 
§ 15.11 We inject fake sources with flux densities ranging from 
—5 to 20 mJy into a synthesised noiseless sky at random 



positions. We iterate this process 10000 times and the prior 
distribution is given as the averaged flux distribution of the 
sources. The deboosted flux densities are given in Table [2] 



5 1.1 MM NUMBER COUNTS 

5.1 Number Counts of the ADF-S 

We create number counts for the 50% coverage region, where 
the noise distribution is more uniform, the survey complete- 
ness at faint flux densities is high, and the number of false 
detections are low compared to the 30% coverage region. 
We employ the Bayesian method, which is now commonly 
used for derivi ng number counts in millimetre / submillimetre 
surveys (e.g.. ICoppin et alJ I2005L I200& IPerera et al.ll2008l : 



lAustermann et alj l2009, 2010l: IScott et al.ll2Q10h 

The PFD of each source is calculated in the same man- 
ner as used in flux deboosting (§ 14.6)) . We adopt the best-fit 
Schecht er function of the SC UBA/SHADES 850 fim number 
counts (jCoppin et all 2006) scaled to 1.1 mm as an initial 
prior distribution function. We create 20000 sample cata- 
logs by bootstrapping off (i.e., sampling with replacement 
the PFDs. We sample only from the PFDs with a crite - 
ria of P(S < 0) ^ 0.05 as adopted in IScott et ail (|2010L 
where P(S < 0) is the probability that the flux densities 
of the detected sources deboosted to <0 mJy. The number 
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Table 3. Differential and cumulative number counts in the ADF- 
S. The flux bin centers for the differential counts (first column) are 
weighted by the assumed prior. The errors are 68.3% confidence 
intervals. 



Table 4. Differential and cumulative number counts in the SXDF 
and the SSA 22 field. The flux bin centers for the differential 
counts (first column) are weighted by the prior. The errors are 
68.3% confidence intervals. 



Flux Density 


Differential 


Flux Density 


Cumulative 


(mJy) 


(mJy- 1 deg- 2 ) 


(mJy) 


(deg- 2 ) 


1.39 
2.41 
3.42 
4.42 
5.43 
6.43 


1074+ 123 
1L " -137 

3121- 
109t£ 
631*1 


1.0 
2.0 
3.0 
4.0 
5.0 
6.0 


15901^ 
204+Ji 

yo -26 
oi+ll 
dl -15 
5 8+ 3 - 4 



of sources in each sample catalog is chosen randomly as a 
Poisson deviate from the real number of sources. The sur- 
vey completeness is calculated by tracing output to input 
sources from the simulated maps. After correcting for the 
survey completeness, the mean counts and the 68.3% confi- 
dence intervals in each flux bin with a bin size of 1 mJy are 
calculated from the 20000 sample catalogs. The deboosted 
fluxes in the source catalogue range from 1 to 6.4 mJy. Each 
of the derived number counts for the 20000 sample catalogs 
is fitted to a Schechter function of the form 



dN 
~dS 



3 mJy 



exp 



-(5 — 3 mJy) 



(i) 



where A^3mj y is a differential counts at 3 mJy and best-fit 
parameters are obtained in S' — iNfemjy parameter space. We 
adopt the above Schechter functional form with a — — 2 
since it well describes nu mber counts derived in previous 
deep SMG surveys (e.g.. ICoppin et al . 2006; IPerera et al.l 
120081 : lAustermann et alJl2009L l2010l : IScott et al]l20ld ). The 
derived best-fit function is then used as a new prior distri- 
bution function and the procedure described above is re- 
peated. The resultant differential and cumulative number 
counts are presented in Fig. Fig. and Table [3] The 
errors indicate the 68.3% confidence intervals. The best-fit 
parameters of the Schechter functional form in equation JT]) 
are iY 3mJy = 169 ± 19 and S' = 1.48 ± 0.20. 



5.2 Number Counts of the SXDF and the SSA 22 
Fields Surveyed by AzTEC/ASTE 

We extract number counts for two other deep fields surveyed 
by AzTEC on the ASTE: the Subaru/XMM Newton Deep 
Field (SXDF) and the SSA 22 field. 

The SXDF is a blank field with deep multi-wavelength 
observations from X-ray to radio. The AzTEC/ASTE obser- 
vations covered ^0.27 deg 2 of the central part of the SXDF 
with an rms noise level of ^0.5-0.9 mJy (S. Ikarashi et al. 
in preparation), which is about a factor of two deeper than 
the A zTEC/JCMT survey of this field (|Austermann et all 
l2010h . In total ~200 sources (^3.5cr) are detected. 

The SSA 22 field is thought to be a proto-cluster region 
since it has an overdensity of UV/optically selected galax- 
ies such as Lyman-Q! emitters ( LAEs) and Lyman - break 
galaxies (LBGs) at z ~ 3.1 (e.g., ISteidel et alJll998l l2000l : 
lHavashino et al. 1 12004 iMatsuda etHI l2005h . The SSA 22 



Flux Density 
(mJy) 



Differential 
(mJy -1 deg -2 ) 



Flux Density 
(mJy) 



Cumulative 
(deg" 2 ) 



SXDF 



1.38 


1110+ 148 

±±±U -148 


1.0 


15781^ 


2.40 


314 tf s 


2.0 


468 +f° 


3.41 




3.0 


154 tit 


4.41 


6Z -14 


4.0 


59 +13 

° y -19 


5.42 


10 +5.7 
10 -9.0 


5.0 


97 +7.6 
Z< -13 


6.42 


74 +3.5 
-6.9 


6.0 


14 +5-1 

L * -8.9 



SSA 22 



1.40 


^79+152 
°' z -178 


1.0 


1006^!^ 


2.42 


9CT1+52 


2.0 


404+59 


3.43 


103+^ 


3.0 


182+g 


4.44 


46 +13 

^° -17 


4.0 


80 + 16 
ou _ 21 


5.44 


19 + 7 - 7 
-li 


5.0 


04 +8.8 
°^ -13 


6.44 


7 q+3.3 
'• U -6.6 


6.0 


15 + 4 " 2 
iD -7.0 



Table 5. Best-fit parameters of parametric fits to differential 
number counts in the ADF-S, SXDF, and SSA 22 fields using 
the Schechter functional form in equation ( [ T ]) . The errors are la 
uncertainty. 



Field 


A^3mJy 


S' 




(mJy -1 deg -2 ) 


(mJy) 


ADF-S 


169 ± 19 


1.48 ±0.20 


SXDF 


144 ± 17 


1.24 ±0.15 


SSA 22 


132 ± 18 


1.85 ±0.38 



field characterized by the overdensity is a good compar- 
ison field with other blank fields to see the relation be- 
tween SMGs and other galaxy populations and the large- 
scale structure of the universe. The AzTEC/ASTE observa- 
tions covered ^0.28 deg 2 with an rms noise level of ^0.6- 
1.2 m Jy, and detected ^100 sources (^3.5cr) (iTamura et al.l 
2009, Tamura et al. in preparation). 

The procedure and parameters in data reductions and 
extracting number counts of these two fields are the same 
as used in this paper. The differential and cumulative num- 
ber counts are presented in Fig. [71 Fig. [HI and Tabled The 
best-fit parameters of the Schechter functional form in equa- 
tion (pQ) are shown in Table [5] 



5.3 Comparison among 1-mm Surveys 

In Figs and [H we compare the number counts in 
the ADF-S, the SXDF, and the SSA 22 fields surveyed 
by AzTEC/ASTE with those of previous 1 -mm survey; 
The AzTEC/JC MT surveys of GOODS- N (|Perera et al.1 
l20Q8h . COSMOS (jAustermann et alJ l2009h , and SHADES 
(combined counts in the Lockman Hole and SXDF) 
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Figure 7. Differential number counts in the ADF-S, the SXDF, 
and the SSA 22 fields surveyed by AzTEC/ASTE. The counts 
are calculated for the 50% coverage regions. The error bars rep- 
resent the 68.3% confidence intervals. The red curve represents 
the best-fit Schechter functional form of the ADF-S: dN/dS = 
N 3 mJy(S/3 mJy)- 1 exp(-(^-3 mJy)/S") with N 3mJy = 169±19 
and S f = 1.48 ±0.20. The number counts of previous 1.1 mm sur- 
yeys are also shown: GO ODS-N dPerera et al.1 120081). COSMO S 
(lAustermann et al.l l2009h. SHAD ES (lAustermann et all l2010h . 
and GOODS-S (jScott et al.l l2010h . The horizontal dashed line 
represents the survey limit of the ADF-S, which Poisson devi- 
ations to zero sources per bin 32.7% . The bin centers of some 
number counts are shifted by ±3% for easier comparison. 



(|Austermann et all lioioh. the AzTEC/ASTE survey of 
GOODS-S (jScott et alfeoioh . and 1.2 mm MAMBO surveys 
of the Lockman Hole and ELAIS N2 (jGreve et al.ll2Q04h . 

The ADF-S and the SXDF provide the tightest con- 
straints on the faint end of the number counts because 
of their depth and large survey areas. On the whole, the 
1 mm counts of various surveys are consistent within errors. 
This is interesting since the SSA 22 field has overdensity of 
UV/optically-selected galaxies. It is possible that the over- 
density of sources at z = 3.1 traced by the UV/optical galax- 
ies does not significantly change the SMG number counts 
given the large volume and redshift-space sampled by the 
mm- wavelength observations. Compared to the ADF-S, the 
counts from GOODS-N and COSMOS are higher, while the 
counts from SHADES are lower. The overdensity of bright 
SMGs in AzTEC/COSMOS compared to other blank fields 
has been shown to be correlated w ith foreground structure 
at z < 1 (|Austermann et al.ll2009h . Since the GOODS and 
the SHADES fields have no known biases, the diversity in 
the number counts likely arises from cosmic variance given 
the small areas of these surveys. 

5.4 Comparison with Models 

We compare the cumulat ive number counts with the m odels 
of iTakeuchi et all (2001a, b). IFranceschini et al.1 (|2010h . and 



Figure 8. Cumulative number counts in the ADF-S, the 
SXDF, and the SSA 22 fields surveyed by AzTEC/ASTE. 
The counts are calculated for the 50% coverage regions. 
The error bars indicate 68.3% confidence intervals. The 
red curve represents the best-fit Schechter functional form 
of the ADF-S. Number counts of previous 1 mm surveys 
are also shown: t he AzTEC/ JCMT surveys of GOO DS-N 
dPerera et al.l l2008h. COSMOS (lAustermann eFaXl I2009T ). and 



S HADES jAustermann et al.ll201oh . the AzTEC/ASTE GOODS- 
S (jScott et al.ll2010h . 1.2 mm MAMBO surveys of the Lockman 
Hole and ELAIS N2 regions flGreve et al.ll2004t ). The horizontal 
dashed line represents the survey limit of the ADF-S, which Pois- 
son deviations to zero sources per bin 32.7%. The bin centers of 
some number counts are shifted by ±3% for easier comparison. 



iRowan- Robinson! |2 



which successfully reproduce the 
observed CIB and number c ounts at IR and submillim etre 
wavelengths. The mode ls of IFranceschini et alJ ((2010) and 
I Rowan- Robinson! (2009) are constructed to match observed 
1.1 mm number counts o f AzTEC surveys of C OSMOS field 
(lAustermann et al]l2009h and GOODS-N field (jPerera et al.1 

120081 ). respectively; 

The model of Takeuchi et al consists of three 

components: (i) the FIR spectral energy distribution (SED) 
based on the IRAS colour- luminosity relation at 60 /xm and 
100 jLtm; (ii) the local 60 /im luminosity function adopted 
from the IRAS data; and (iii) galaxy evolution with redshift. 
The 60-/xm luminosity of a galaxy is thus described as a 
function of redshift, assuming pure luminosity evolution: 



L 60 (z) = L 60 (z = 0)f(z), 



(2) 



where Leo is the luminosity at 60 /im and Lqo(z = 0) rep- 
resents the local 60 fim luminosity functio n. The form of 
f(z) is a stepwise nonparametric function. I Takeuchi et al.l 
( 200 laf) assume three evolutionary scenarios (see Figure 2 of 
ITakeuchi et al.1 l2001ah within the permitted range derived 
from the observed CIB and number counts at 15, 60, 90, 
170, 450, and 850 /im: (i) Evolution 1: f(z) rises steeply 
between z — — 0.5, reaching f(z) = 10, is constant from 
z = 0.5 — 2, and decreases slowly between z = 2 — 5; (ii) 
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Evolution 2: f(z) quickly rises between z = — 0.5, peaking 
with f(z) = 30 between z = 0.5 — 0.75, and having a long 
plateau with f(z) = 10 between z = 0.75 — 5; (iii) Evolu- 
tion 3: f(z) rises between z = — 1, peaks with f(z) = 30 
between z — 1 — 2, and has a plateau with f(z) = 20 between 
z = 2 — 5. They adopt two additional evolution models which 
are made by modifying Evolution 1 to rise to f(z) = 10 1-5 
and f(z) = 10 2 between = 1-2, which we will refer to as 
"Evolution 4" and "Evolution 5" , respectively. Figure[9]com- 
pares 1.1 mm observed number counts to the five models, 
along with a no-evolution model. The no-evolution, Evolu- 
tion 1, and Evolution 2 models are lower than the observed 
counts, suggesting that significant luminosity evolution with 
f(z) > 10 is needed. Given that these five models fail to ex- 
plain the observed counts, the models need to be modified. 
Possible ways to solve this discrepancy are (i) to change 
the functional form of the evolutionary scenarios, and (ii) 
to make the luminosity evolution dependent on luminosity. 
These issues are discussed in more detail in Takeuchi et al., 
in preparation. 

The model of iFranceschini et al.l (|2010h is constructed 
to reproduce the latest observed counts at 15, 24, 70, 350, 
850, and 1100 /im, the redshift dependent luminosity func- 
tions at 15 /im, and the CIB. They assume both lumi- 
nosity and number density evolution, and create number 
counts starting from the IRAS 12 /im luminosity function. 
The model population consists of four galaxy classes: non- 
evolving normal spirals, type-I AGNs, starburst galaxies 
of moderate luminosities (or LIRGs), and very luminous 
starburst galaxies (or ULIRGs). The four populations fol- 
low different evolution in lumino s ity an d number density. 
The model of IFranceschini et al.l (|2010h overestimates the 
ADF-S counts, while it is consistent with the counts of 
the AzTEC/COSMOS survey. This is not surprising, since 
the model is created to reproduce the 1.1 mm counts of 
the AzTEC/COSMOS, wh ere a significant overd ensity of 
sources has been re ported ( Austerm ann et al1l2009l ). 

The model of IRowan- Robinsonl (200 9|) is created by 
modifying the model of IRowan- Robinsonl feOOll ) to repro- 
duce the latest observed counts, particularly at 24 /im. 
The model assumes pure luminosity evolution. The model 
consists of four spectral components: infrared cirrus, M82- 
like starburst, Arp 220- like starburst, and AGN dust torus. 
I Rowan- Robinsonl (2009) creates three models with forma- 
tion redshift s of Zf = 4, 5 and 10. Figure [TUl shows that the 
model with Zf = 4 well describes the ADF-S counts down 
to the faint end, but at the bright end it overestimates the 
ADF-S counts. 

All of the models presented in this section do not match 
simultaneously the faint and bright end of observed counts, 
requiring the models to be modified. 



6 CONTRIBUTION TO COSMIC INFRARED 
BACKGROUND 

We estimate the fraction of the CIB resolved by the ADF- 
S survey. The total deboosted flux density of the ^3.5cr 
sources in the 50% coverage region is 1.9 =b 0.03 Jy deg -2 . 
The expected 1.1 mm background as measured by the 
Cosmic Background Explorer satell ite is 18-24 Jy deg -2 
(|Puget et al.ll 19961 : iFixsen et al.ll 19981 ) , therefore we have re- 




Flux Density S (mJy) 



Figure 9. C omparison of 1. 1 mm o bserved number counts with 
the models of lTakeuchi et all (|2001al V Descriptions of the models 
are in § 15.41 




1 



Flux Density S (mJy) 

Figure 10. Co mparison of 1.1 mm obse rved number counts with 
the m odels of Rowan- Robinson (2009) and IFranceschini et all 
(2010). Descriptions of the models are in § 15.41 



solved about 7-10% of the CIB into discrete sources. This is 



similar to the resolved fraction of the CIB re 


ported by other 


1 mm blank field survevs (IGreve et al. 2004 


; 


Laurent et al.l 


2005: iMalonev et al.1 120051: IScott et al.l 2008 


. l2010h. This 



could be caused by the survey incompleteness due to the 
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confusion noise and the fewer bright sources in the ADF-S, 
despite the deeper sensitivity compared to the other surveys. 

To estimate the total integrated flux density corrected 
for the survey incompleteness, and to include fainter sources 
below the detection threshold, we integrate the best-fit 
Schechter function of the ADF-S obtained in § [5] The inte- 
gration of the best-fit function at ^1 mJy, where the num- 
ber counts are tightly constrained, is 2.9 Jy deg -2 , which 
corresponds to 12-16% of the CIB at 1.1 mm. This suggests 
that a large fraction of the CIB originates from submm-faint 
sources for which the number counts have not yet been con- 
strained. Integration of the best-fit number counts extrapo- 
lating to lower fluxes (down to mJy) results in a total flux 
density of 5.7 Jy deg" 2 , which is only 24-32% of the CIB 
at 1.1 mm, suggesting that the faint-end slope of the actual 
number counts should be steeper than that of the present 
best-fit model. It is possible that a Schechter functional form 
is not appropriate for representing 1.1 mm number counts. 



7 REDSHIFT CONSTRAINT 

In order to constrain redshifts of AzTEC sources, we com- 
pare AzTEC sources in the 30% coverage region with the 
far-infrared images obtained by AKARI/FIS (Shirahata et 
al. in prep.). The ADF-S is observed at four bands: 65, 90, 
140, and 160 /mi, with FWHMs of 37", 39", 58", and 61", 
respectively. The detection limits of the AKARI data are 
46.5, 15.7, 183, 608 mJy (3a) at 65, 90, 140, and 160 /im, re- 
spectively. We compare the AzTEC sources with the 90 /xm 
source catalog, which is most sensitive and reliable among 
the four bands, and found only 11 AzTEC sources are within 
a 20" radius from the 90 /im sources. A detailed multiwave- 
length study of these sources will be made in a future paper. 
We constrain the redshifts of the AzTEC sources using their 
flux ratios between 1.1 mm and 90 /im. Figure [TT] shows the 
expected flux ratio as a function of redshi ft for two different 
SED models: Arp 220 (|Silva et al.lll998T ). a typical ultra- 
luminous infrared galaxy (ULIRG), and the average SEP 
of 76 SMGs with spectroscopic redshifts (Mic halowski et aT] 
l2010h . The horizontal dotted lines represent the flux ratios 
of AzTEC sources with 90 /im counterpart candidates. The 
shaded region represents upper limits on the flux ratios of 
the remaining AzTEC sources from the 3a detection limit at 
90 /xm. The figure suggests that most of the AzTEC sources 
are likely to be at z > 1.5. 

At z > 1, the flux density of SMGs is nearly redshift 
independent and thus is proportional to the IR luminosity. 
By scaling the IR luminosities of the SED models, we es- 
timate IR luminosities of the AzTEC sources to be ^3-14 
xlO 12 L©. If the emission is powered solely by star forma- 
tion activity, th e inferred SFRs are -500-2400 M yr" 1 
(jKennicuttlll998T ). 



8 COSMIC STAR FORMATION HISTORY 
TRACED BY 1.1 MM SOURCES 

8.1 Star- format ion Rate Density 

From UV/optical surveys, the total SFR per unit comov- 
ing volume (SFR density) is observed to increase with 
redshift from z ~ 7 to z ~ 3, peak at z ~ 3-1, and 




Figure 11. Observed flux ratios between 1.1 mm and 90 fim. 
as a function of redshift. The solid and dashed curves represent 
the SED models of Arp 220 (|Silva et al.lll998 | ) and the average 
SED of 76 SMGs with spectroscopic redshifts I Michalowski et al.l 
2010), respectively. The horizontal dotted lines indicate the flux 
ratios of AzTEC sources with 90 /im counterpart candidates. The 
shaded region represents the 3a upper limits of 90 fim for AzTEC 
sources without 90 fim detections. 



decline steadily from z — 1-0 (e - g.. iMadau et al.. 
i Lillv et al.lll996l: ISteidel et al.l [l999l : iGiavalisco et al.l 



1996 



2004 



Bouwens et al. uOloh . It is suggested that SFRs per unit co- 



moving volume (SFR density) increases with redshift from 
z ~ 7 to z ~ 3, peaks at z ~ 1-3, and decreases from 
z ~ 1. However, the SFR density derived by UV/optically 
selected galaxies has large uncertainty due to the extinction 
by dust. It is also possible that dusty galaxies are missed 
entirely by these previous studies. Comparatively speaking, 
millimetre and submillimetre wavelengths have a great ad- 
vantage in tracing dusty starburst galaxies at high redshifts. 
Previous submillimetre surveys suggested that SMGs con- 
tribute significan tly (> 10-20%) to the cosmic SFR densit; 

at z ~ 2-3 (e.g., iHu ghes eT^ilT" 

lAretxaga et al.ll2007l : iDve et al.ll2008l : IWardlow et al.ll2010h 
We estimate the SFR density contributed by 1.1 mm 
sources using the best-fit number counts in the ADF-S de- 
rived in §[5] FIR luminosities of the 1.1 mm sour ces are cal- 
culated by assuming the SED models of Arp 220 (ISilva et al 



ty 

'5; 



1998) and the average SED of SMGs dMichalowski et al 



l2010f ), and SFRs are derived from FIR luminosity using the 
equation of lKennicuttl (|l998h . The largest uncertainty comes 
from the lack of redshift information. Since the redshifts of 
the 1.1 mm sources are not known, we assume redshift distri- 
butions based on previous studi es. The largest s pectro scopic 
sample of SMGs obtained bv IChapman et ahl (|2005l ) has a 
median redshift of z = 2.2 with an interquartile range of 
1.7-2 .8, and the r edshi ft distribution is well fitted by Gaus- 
sian. IPope et aL (2006) found a median redshift of z — 2.2 
with an interquartile range of 1.4-2.6 using spectroscopic 
and photom etric redshifts o f SMG s in the Hubble Deep 
Field-North. lAretxaga et al.l ([20071 ) estimated photometric 
redshifts of SHADES sources and found a median of z — 2.4 
with an interquartile range of z = 1.8- 3.1, and that redshift 
distribution has a near- Gaussian form. IChapin et all <|2009h 
derive a higher median redshift of z = 2.7 using spectro- 
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Table 6. Comoving star formation rate density averaged in red- 
shift bins estimated from the best-fit number counts in the ADF- 
S. Four redshift distributions are assumed. The SFR densities are 
calculated by integrating the number counts at ^1 mJy and using 
the SED model of Arp 220. 

(zcjCz) Comoving Star-formation Rate Density 

(M© yr _1 Mpc -3 ) 





z = 1-2 


z = 2-3 


z = 3-4 


z = 4-5 


(2.2, 0.5) 


3.0 x 10~ 2 


4.6 x 10~ 2 


4.1 x 10~ 3 


1.2 x 10~ 5 


(2.4, 0.5) 


1.9 x 10~ 2 


5.1 x 10- 2 


8.5 x 10- 3 


5.2 x 10- 5 


(2.4, 1.0) 


2.3 x 10- 2 


2.9 x 10- 2 


1.6 x 10- 2 


3.8 x 10- 3 


(2.7, 1.0) 


1.7 x 10- 2 


2.9 x 10- 2 


2.1 x 10- 2 


6.5 x 10-3 




scopic and photometric redshifts of 1.1 mm sources detected 
in the AzTEC/JCMT GOODS-N survey. 

Based on these measurements, we adopt Gaussian func- 
tional forms of redshift distributions with central redshifts of 
z c = 2.2, 2.4, and 2.7. We assume a z — 0.5 and a z — 1.0 for 
narrow and broad redshift distributions, respectively. These 
redshift distributions are consistent with the fact that the 
most of the AzTEC sources are at z > 1.5 (see § [7J). We 
calculate the total flux by integrating the differential num- 
ber counts at ^1 mJy and distribute the total flux following 
the assumed redshift distributions. If we integrate the num- 
ber counts down to 0.1 mJy, the total flux density would 
increase by about a factor of two. 

The derived SFR densities using the Arp 220 SED 
model are shown in Figure [12] and the average values in 
redshift bins are presented in Table \6\ The results de- 
rived from the two assumed SED models are consistent 
within 30%. Although the derived SFR densities largely de- 
pend on the assumed redshift distribution, they are within 
the range of those derived in previous studies with SMGs 
(iHughes et al J 1 19981: [Chapman et al.1 20051; IXretxaga et"ai] 
120071 : iDve et al.ll2008l : IWardlow et alJl201Ch . In Figure [12] 
extinction-corrected SFR densities derived from previous 
UV/optical o bservations are also shown for comparison 
(|HoDkinj|2003 ). The SFR densities of 1.1 mm sources are 
lower by about a factor of 5-10 at z ~ 2-3 compared to that 
of the UV/optically-selected galaxies. 

In § 1 we found that integrating the AzTEC/ADF-S 
number counts at ^1 mJy accounts for 12-16% of CIB at 
1.1 mm. If we assume that the rest of the CIB comes from 
fainter sources, the total SFR density contributed by 1.1 mm 
sources including <1 mJy sources would increase by about 
a factor of 6-8, which is comparable to or higher than that 
of the UV/optically-selected galaxies at z ~ 2-3. We note 
that in this case the faint 1.1 mm sources and UV/optical 
sources can overlap. The large contri bution of dusty g alaxies 
to the SFR density is suggested bv lGoto et al.l (|2010h based 
on 8 /im and 12 /mi observations. They found that the dust- 
obscured SFR density at z ~ 2 is -0.5 M yr" 1 Mpc -3 , 
which is consistent with the SFR density of 1.1 mm sources 
including the fainter (<1 mJy) population sources. 



8.2 Stellar Mass Density 

We estimate the fraction of the stellar mass in the present- 
day universe produced by 1.1 mm sources with ^1 mJy by 



Figure 12. Star-formation rate density estimated from the best- 
fit number counts in the ADF-S. Four redshift distributions are 
assumed. The SFR densities are calculated by integrating the 
number counts at ^1 mJy and using the SED model of Arp 220. 
The results from previous UV-optical obser vations are pl otted 
for comparison (compiled data are taken from Hopkins 2004, and 
references therein). 



integrating the SFR density derived in the previous section. 
The present-day stellar m ass density is estimated from local 
luminosity functions fe.g.JCole et al.ll200ll : iBell et alJliooi : 
iKaiisawa et al.|[2009h . ICole et al.l (|200lh derived the present- 
day stella r mass density o f (5.6 ± 0.8) xlO 8 M Mpc -3 as- 
suming a ISalpeterl (|1955h initial mass function. The time 
integration of the SFR densities at z > 1 for the four 
assumed redshift distributions ((z c ,cr z ) =(2.2, 0.5), (2.4, 
0.5), (2.4, 1.0), and (2.7, 1.0)) yields 1.2 x 10 8 M Mpc" 3 , 
1.0 x 10 8 M Mpc" 3 , 1.0 x 10 8 M Mpc -3 , and 0.90 x 
10 8 M Mpc -3 , respectively. This corresponds to -20% of 
the present-day stellar mass density. This is an upper limit 
since materials forming massive stars returned to the inter- 
stellar medium via stellar winds and supernovae explosions. 
The fraction of stellar mass returned to the ISM, called as 
the r ecycled fraction, i s estimated in sem i - analytical models 
(e.g.JCole et alll2000l : baueh et al.ll2005l : lLacev et al.ll2010l : 
iGonzalez et al llpOlOh . and it depends on the IMF: 0.41 for 
thelKennicutti (Il983l) IMF an d 0.91 for a top-heavy IMF 
ICole et all (|2000h : lLacev et al.l (|2010h . If we assume the re- 
cycled fraction of 0.41 and 0.91, the fraction of the stel- 
lar mass in the present-day universe produced by 1.1 mm 
sources with ^1 mJy decrease to —10% and a few percent, 
respectively. 



9 SUMMARY 

We performed a 1.1 mm deep survey of the AKARI Deep 
Field South (ADF-S) with AzTEC mounted on the ASTE, 
obtaining one of the deepest and widest maps at millimetre 
wavelengths. The 30% and 50% coverage regions have areas 
of 909 and 709 arcmin 2 , and noise levels of 0.32-0.71 mJy 
and 0.32-0.55 mJy, respectively. We detected 198 previously 
unknown millimetre-bright sources with 3.5-15.6cr in the 
30% coverage region, providing the largest 1.1 mm source 
catalog from a contiguous region. 
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We constructed differential and cumulative number 
counts in the ADF-S, the SXDF, and the SSA 22 field which 
probe fainter flux densities (down to 1 mJy) compared to 
previous surveys except for the AzTEC/ASTE GOODS-S 
survey. On the whole, the 1 mm counts of various surveys 
are consistent within errors. We compare the number counts 
with t he luminosity evolut i on m odels ofjTakeuchi et al. 
(2001a|), iFranceschini et all (l2010h . and I Rowan- Robinson! 
( 20091 ) . Comparison with the iTakeuchi et al. (2001a) model 
suggests that a luminosity evolution with a factor >10 is 
needed to explain the observed num ber counts. The observed 
number counts favor the model of Rowan- Robinson (2009) 
with Zf = 4, but none of these models simultaneously match 
both the bright and faint end of the number counts from 1— 
10 mJy. 

In the ADF-S survey, we resolve about 7-10% of the 
CIB at 1.1 mm into discrete sources. The integration of 
the best-fit number counts in the ADF-S down to 1 mJy 
reaches 12-16% of the CIB. This suggests that the large 
fraction of the CIB at 1.1 mm originates from faint sources 
(ft.imm < 1 mJy) for which the number counts have not 
yet been constrained. The integration of the best-fit number 
counts extrapolating to mJy accounts for only 24-32% of 
the CIB, suggesting that the faint-end slope of the number 
counts is steeper than that given by our best-fit model. 

The redshifts of the AzTEC sources are constrained 
from their flux ratios between 1.1 mm and 90 fim from the 
AKARI/FIS. Most of the AzTEC sources are not detected 
at 90 /im, suggesting that they are likely to be at z > 1.5. 
Assuming z > 1, the inferred IR luminosities of the AzTEC 
sources are -(3-14) xlO 12 L©, and their SFRs inferred from 
the IR luminosities are -500-2400 M yr -1 . 

We derived the cosmic SFR density contributed by 
1.1 mm sources using the best-fit model to the differential 
number counts. Although the derived SFR density largely 
depends on the assumed redshift distribution, our estimates 
are within the range of those derived in previous studies with 
SMGs. The SFR density of 1.1 mm sources with ^1 mJy at 
z — 2-3 is lower by about a factor of 5-10 than those of 
UV/optically-selected galaxies. If we consider the fact that 
the contribution of 1.1 mm sources with ^1 mJy to the CIB 
at 1.1 mm is 12-16%, the SFR density of 1.1 mm sources 
including those fainter than Si.imm < 1 mJy would become 
comparable to or higher than those of UV/optically-selected 
galaxies. The fraction of the present-day stellar mass of the 
universe produced by 1.1 mm sources with ^1 mJy at z ^ 1 
is —20%, calculated by the time integration of the SFR den- 
sity. If we consider the recycled fraction of 0.41 and 0.91, the 
fraction of stellar mass produced by 1.1 mm sources becomes 
— 10% and a few percent, respectively. 
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